arXiv:1506.08763v2 [quant-ph] 22 Sep 2015 


Quantum Zeno effect in parameter estimation 

Alexander Holm KiilericlQ and Klaus M0lmer 
Department of Physics and Astronomy, Aarhus University, 

Ny Munkegade 120, DK 8000 Aarhus C. Denmark 
(Dated: September 23, 2015) 

The quantum Zeno effect freezes the evolution of a quantum system subject to frequent measure¬ 
ments. We apply a Fisher information analysis to show that because of this effect, a closed quantum 
system should be probed as rarely as possible while a dissipative quantum systems should be probed 
at specifically determined intervals to yield the optimal estimation of parameters governing the sys¬ 
tem dynamics. With a Bayesian analysis we show that a few frequent measurements are needed to 
identify the parameter region within which the Fisher information analysis applies. 
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I. INTRODUCTION 

Quantum systems, ranging from atomic systems to 
field modes and mechanical devices are useful precision 
probes for a variety of physical properties and phenomena 
DU]. When we deal with a single quantum system, we 
must take into account that the measurements by which 
we extract information about its evolution yield random 
results and that they cause a measurement back action 
on the system. For some problems this back action may 
be favorable as it randomly quenches the system and thus 
triggers a transient evolution with temporal signal cor¬ 
relations and these may depend more strongly than the 
steady state on the desired physical properties Dili. 

For too strong or too frequent measurements, however, 
the back action may completely dominate the evolution. 
This is manifested in the Quantum Zeno effect (QZE) 
D, named after Zeno’s arrow paradox, which inhibits 
population transfer between discrete states in a quantum 
system due to frequent observations. The effect was first 
illustrated in ion trap experiments [3] and has since been 
used to account for experiments on many other systems, 
e.g DlHj. 

The intuition behind the quantum Zeno effect has stim¬ 
ulated proposals for application in quantum information 
processing |9]i 6.g. for entanglement protection |in| . 
preservation of systems in decoherence free subspaces 
Dim and error suppression in quantum computing |12| . 
In many of these studies the measurement back action 
is actually not unambiguously identified as the crucial 
element, but experiments and calculations have clearly 
established the intuition behind the quantum Zeno effect 
as a useful and inspirational source for new ideas and 
proposals. 

In this study, we focus on information retrieved di¬ 
rectly by the measurements and, closer to the original 
Zeno paradox, we ask to what extent frequent measure¬ 
ments on a quantum system prevent the observation - 
through the same measurements - of its dynamical evo¬ 
lution. This analysis is pertinent for the use of quantum 
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systems as sensitive probes, and we shall quantify the 
information content in measurement records and derive 
theoretical expressions for the Fisher information m , to 
reveal the extent to which high precision probing is in¬ 
hibited by the quantum Zeno effect in the limit of strong 
frequent measurements. The Fisher information quanti¬ 
fies the asymptotic precision in any estimate and has in 
recent works been applied to analyze quantum parameter 
estimation in various schemes of repetitive and continu¬ 
ous measurements, see e.g. DanMii. 

As our main example, we consider the estimation of 
the Rabi frequency SI of a driven two-level system sub¬ 
ject to projective eigenstate measurements. A similar 
setup is considered in ref. m- Looking at FIG. it 
is intuitively clear that the Rabi oscillation dynamics is 
prevented as the interval r between measurements de¬ 
creases, and if we measure too often, the signal holds no 
information about SI. If, on the other hand, we measures 
too rarely we may obtain very little data in any finite 
probing time T. So, what is the optimal value of t for 
the purpose of parameter estimation? Our expression for 
the Fisher information provides the t maximizing the in¬ 
formation in any record obtained in a finite time, and by 
employing a Bayesian inference protocol |18| to the mea¬ 
surement outcomes we discuss further the convergence 
of parameter estimation based on such optimal probing 
intervals. 

The article is organized as follows. In Sec. |TT| we 
present a statistical analysis of the Fisher information 
based on records of projective measurements performed 
at equidistant times. In Sec. |III| we give a simple ar¬ 
gument for the QZE in open systems. In Sec. EYi the 
QZE and Fisher information are discussed in the case 
of two-level systems. Sec. jV] exemplifies the theory in 
terms of a two-level model with and without inherent de¬ 
phasing. In Sec. I VI I a Bayesian parameter estimation 
protocol is applied and strategies that reach the asymp¬ 
totic resolution offered by the Fisher information analysis 
are investigated. In Sec. |VII[ we conclude and discuss 
the relevance of the results for different probing schemes. 



II. PROJECTIVE MEASUREMENTS AND 
FISHER INFORMATION 
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The unconditional evolution of an open quantum sys¬ 
tem is in the Born-Markov approximation given by a 
master equation dtp{t) = Cp(t), with a Liouvillian su¬ 
peroperator on Lindblad form, 

C[p] = -i[H,p] + J2 UkPci - ^{ci^k,p}j , (1) 

k ^ ^ 

where h = 1, and Ck represent different relaxation pro¬ 
cesses. For time independent operators the solution is 
formally given by exponentiation, 

p(t) = e^*[p(0)]. (2) 

If one monitors the time evolution by frequent (contin¬ 
uous) observation the back action conditions the system 
evolution on the measurement outcomes and the evolu¬ 
tion becomes stochastic. 

A realistic approach to continuous or frequent mea¬ 
surements should take the finite bandwidth and noise 
properties of the measurement device into acount, see, 
e.g., HH, but we shall restrict our attention to the ideal 
case of instantaneous, accurate measurements, repeated 
at regular intervals t. The continuous regime approxi¬ 
mated by the limit t —0 thus assumes a corresponding 
high bandwidth of detection. 

While measurement theory can be formulated more 
generally, for simplicity we restrict our analysis to projec¬ 
tive measurements. Each measurement outcome is thus 
an eigenvalue A of a given operator A, occurring with the 
probability 

P(A)=Tr(nAp(f))=mW, (3) 


where !!> = |A) (A| is a projector on the corresponding 
eigenstate. The projection postulate states that the con¬ 
ditional state after a measurement performed on the state 
pit) at time t is |2[)| . 


PA(t) 


nA/o(0riA 

P(A) 


(4) 


A. Fisher information of projective measurement 
records 


The acquisition of information about unknown physi¬ 
cal parameters from a measurement record follows Bayes’ 
rule: Given the stochastic measurement outcome D, the 
probability that an unknown parameter has a given value 
0, is given by the corresponding outcome probability con¬ 
ditioned on the value 9 and its unconditional (prior) prob¬ 
abilities. 
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FIG. 1. (Color online) Illustration showing the excited state 
population as a function of time for a two-level system driven 
on resonance, and randomly projected at regular intervals ac¬ 
cording to measurements in the eigenstate basis. The mea¬ 
surement intervals are different in the different panels, (a): 
T = 2.5f2“^, (b): T = 1.75f2“^, (c): r — 0.75I2~^ and (d): 
T — 0.25f2~^, and show the Quantum Zeno suppression of the 
coherent Rabi oscillations. 


The Cramer Rao bound (CRB) |21| . 

|AS(0)1» > (6) 

yields a lower bound of the statistical variance [A5'(0)]^ 
on any unbiased estimate of the unknown quantity 9 
(with the true value 9o). Here, 

F(9) = -E[d^lnP(Dl9)], (7) 

is the Fisher information which quantifies the dependence 
of the data outcome probabilities on the unknown pa¬ 
rameter. The Fisher information and the Cramer Rao 
bound yield the asymptotic precision of the best possible 
estimate, and normally assumes a large number of repe¬ 
titions JV of the experiment (and a corresponding factor 
1/N in Eq. In our analysis a self-averaging occurs 

due to the repetitive measurements, and we shall find 
the corresponding dependence of P(9) on JV and, equiv¬ 
alently, on the total time of the probing T. 

A data record of V = T/t projective measurement 
outcomes obtained at times tj = jr, D = {Aj}^]^, occurs 
with the probability 

P(D|0) = nTr(|A,)(A,|e^"[p(t,_i)]) 

3 

3 

assuming that the initial state is one of the eigen¬ 
states, p(0) = |Ao)(Ao|. Since by Eq. Q the sys¬ 
tem is always projected (reset) in one of the eigen¬ 
states upon detection, the probability for an entire data 
record factors into a product of conditional probabilities 
P(A,|A,_i,0) = Tr(|A,)(A,|e^"[|A,_i)(A,_J]), where 
we recall that the time evolution operator Q depends 
on 9. The data record is fully represented by the set of 
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numbers nim counting the occurrences of subsequent de¬ 
tection in states (Aj_i,Aj) = The mean num¬ 

ber of such events during a total of N measurements 
is E[nir„] = NP{Xjn\Xi,6)P{\i\9), where the probability 
that the average measurement will yield A; may be calcu¬ 
lated as P{Xi\9) = pf x ’"^hh the stationary (steady 
state) solution of the non-selective evolution in Eq. 
where the average measurement back action is included 
as a dephasing of the atomic coherence at a rate oc 
(See Sec. VB below). 


The conditional probability for subsequent measure¬ 
ments, P{Xj = Am|Aj_i = Xi,9) does not depend on j, 
and with a total number of measurements N = 
the probability in Eqs. for the data record D — 

{nim} is a multinomial distribution, 


r = T/N, the probability that all measurements yield 
the initial state reads 

N 

/TA 

1 


pquadratic 

^0 


{N)c^ 




(13) 


In the limit N —)■ oo, tends to 1: The QZE 

freezes the system dynamics. 

At small but finite r, it is instructive to set (T/7V)^ = 
t{T/N) in Eq. (131 and keep r constant. Then, taking 
the limit qq) = exp(r6T). The proba¬ 

bility to find the system in its initial state thus decays 
on the Zeno time scale r^, given by = —rb. For 
closed system dynamics governed by a Hamiltonian H, a 
vanishes, but we obtain — {H) ) |22| . 


P{D\9)=Y[P{Xm\Xi,9y 


(9) 


Im 


Utilizing J2rn^i^m\Xi,9) = 1, we readily obtain the 
Fisher information for a sequence of N projective mea¬ 
surements, which we can represent as 


FA9) 

N 


= E 


Im 


idePiXm\Xi,9)y 

P{Xm\Xl,9) 


P{Xi\9). (10) 


The result Eq. (101 applies to any parameter estimate 


from a record of projective measurements performed at 
constant intervals r. 

The quantity of relevance in this article is the Fisher in¬ 
formation obtained for a given, long, interrogation time, 
T = Nt, and to address the asymptotic precision as well 
as the optimal value of r, we shall consider the Fisher 
information per time Fr{9)/T — Fr{9)/NT. 


III. ZENO INHIBITED EVOLUTION 

Following Peres |22], we give now a simple account of 
the Quantum Zeno effect. Let po denote a pure state 
of a quantum system subject to evolution by the master 
equation Eq. ([^. The probability that the system will 
be found in its initial state after a short time r is 


Pq — 1 + or -I- hA 


( 11 ) 


where a = Tr(po>C[po]) and b — iTr (po^^[po])- In cases 
where a yf 0, the dominating contribution to Pq is linear 
in r, and the survival probability in the initial state after 
N time steps and projective measurements is 


plinear 


{N)^ 


T 


n N 


( 12 ) 


yielding an exponential decay law Pg'"®®'’'(iV —?> oo) = 

exp (aP) for large N. _ 

In a number of cases, a = 0 (See Sec. |IV| below). Then, 
if the measurement is performed N times at intervals 


IV. QZE AND FISHER INFORMATION FOR A 
PROBED TWO-LEVEL SYSTEM 


We shall now investigate how the Zeno effect af¬ 
fects precision probing by considering a two-level system 
driven by a laser field. In a frame rotating with the fre¬ 
quency of the laser field, the Hamiltonian reads 

H = -5\e){e\ + ^{\e){g\ + \g){e\), (14) 

where 5 is the atom-laser detuning and Q is the laser-Rabi 
frequency. 

Following the previous discussion, projective measure¬ 
ments in the atomic eigenstate basis, Hg = \g) {g\ and 
He = |e) (e|, are performed at intervals r. The condi¬ 
tional, stochastic evolution is exemplified in FIG. [T] for 
different durations of the interval r between measurement 


during which the evolution given by (141 is unitary. De¬ 


creasing T gradually inhibits the coherent dynamics, and 
eventually leads to quantum jumps occurring at the rate 
jA. 

To determine the Fisher information in Eq. ( [Io| , we 
note that the Xj take only two values corresponding to 
measuring the ground (Ay = g) or excited (Ay = e) state 
respectively. In this case. 


FA9) 

N 


P{g\9)(dePig\g,9)A 
Pi9\9A) A - Pi9\9A)) 
P{e\9){deP{e\e,9)y 


+ 


P(e|e,6»)(l-P(e|e,6l))’ 


(15) 


where P{9\9) and P{e\9) are the populations of the 
ground and excited state, in the unconditional steady 
state of the system. Absent energy relaxing processes 
these are both 1. 

Due to the symmetry of the problem, the only relevant 
information in the measurement record concerns how of¬ 
ten consecutive measurements yield identical or different 
results, and we have 


FAQ) _ 

N P{g\g,9){l-P{g\g,9))-' 


{doP{g\g,9)A 


(16) 
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which can be recognized as the Fisher information for 
binomially distributed data. 


For very short durations between measurements 
P{g\g,9) is on the form of Eq. (Ill with po = Iff) (ffl) 
and from Eq. the Fisher information is 


Fr{9) ^ {dga)'^T'^ + 2{dea){deb)T^ + {deb^T^ 
N 


ar + (a^ + b)T‘^ 


(17) 


If a = 0 the QZE freezes out the system dynamics in the 
limit of r —)■ 0, and Fr{9)/T ~ T{dgb)^/b —>■ 0, so despite 
the growing number of measurements, N = T/t ^ oo 
the data holds vanishing information on 9. If, on the 
other hand, a ^ 0, we find a constant Fisher information 
per time in the limit of r —>■ 0, Fr{9)/T ~ {dgaY/a. 
Note, that if a is independent of 9 the Fisher information 
may still vanish even though the dynamics is not frozen 
out for r —>■ 0. 

For simplicity, we restricted this discussion to the case 
of a two-level system, but the relation between a vanish¬ 
ing linear term in the survival probability Eq. (11) and 


Quantum Zeno inhibited parameter estimation equally 
applies to the general case of Eq. (10). 
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V. EXAMPLES 

In this section we exemplify the Zeno inhibited Fisher 
information for Rabi-frequency estimation in a driven 
two-level system. We first discuss the case of unitary 
dynamics between measurements, and we then turn to a 
model including dephasing. 


FIG. 2. (Color online) Upper panel: The Fisher information 
FV(fl) = tT for Rabi frequency estimation in a resonantly 
driven two-level system as a function of the total probing time. 
Results are shown for different durations of the intervals be¬ 
tween measurements. Lower panel: The Fisher information 
Eq. (191 per time for Rabi frequency estimation as a func¬ 
tion of the duration of the intervals between measurements 
assuming r T. Results are shown for different detunings. 


A. Driven system with no dephasing 


For the two-level system with Hamiltonian Eq. (141 the 
time evolution may be solved analytically. 


P(ff|ff,fl) = i 


1-f 


-h cos XT 
^2 


(18) 


where x = + is the generalized Rabi-frequency. 

By Eq. ([^ the corresponding Fisher information per 
measurement for Rabi-frequency estimation is 


Frm 

N 


cos ^ -I- 25“^ sin 


- UF cos^ ^) 


(19) 


In the upper panel of FIG. we show how the Fisher 
information grows as a function of the total interrogation 
time T = Nt for different values of r. Information is only 
obtained when measurement outcomes are acquired, and 
the Fisher information therefore grows in steps. While 
the frequency of the steps increases when the durations 
between measurements are decreased, the figure clearly 


shows a reduction of both the height of the individual 
steps and the effective slope of the Fisher information as 
function of time T. In the lower panel of FIG. the 
slope Fr/T is shown as function of r for different val¬ 
ues of the laser detuning and t T . On resonance 
(S — 0) Eq. (191 reduces to the Fisher information per 
measurement, Ft(UI.)/N = r^, independent of the actual 
value of 0 (equivalent to a linear relationship between 
Fr{Ut)/(TUl~'^) and Or). For a fixed total time, the Rabi 
oscillation dynamics is frozen out and as anticipated by 
the discussion following Eq. 0 any 0-estimate is inhib¬ 
ited by the QZE in the limit r —>■ 0. 

To understand the behaviour for large r, we consider 
the numerator in Eq. (16l. In the upper panel of FIG.[^ 
P(g\g, O) is shown for two slightly different values of O. 
The two candidate values lead to solutions evolving at 
different frequencies acquiring the same phase after a cer¬ 
tain time. This effect contributes the envelope to the 
difference quotient seen in the lower plot (dashed, or¬ 
ange curve). The smaller the difference between the Rabi 
frequencies, the later the peak appears in the envelope 
function, and in the limit of infinitesimally close Rabi 
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FIG. 3. (Color online) Upper panel: The probability of find¬ 
ing a two-level system driven on resonance in the ground 
state as a function of time for Hi = SI (blue curve) 
and SI 2 = 1.1 SI (dashed, red curve). Lower panel: The 
dashed, green curve shows the squared difference quotient, 
[P{g\g, H 2 ) — Pig\g, Hi)] / — SIi], between the probabili¬ 

ties in the upper panel while the orange curve shows the 
squared differential quotient, dnP(g\g, fl), corresponding to 
the limit of infinitesimally close Rabi frequencies. 


frequencies, their Rabi oscillations will only come back 
in phase after an infinitely long time. As the Fisher in¬ 
formation relates to infinitesimally different hypotheses, 
a single measurement at time r = T provides the optimal 
distinction, and the limiting procedure of preparation of 
an initial state followed by a single read out at the final 
time is favoured. Indeed, the quadratic dependence on 
T, causing the Zeno suppression for short measurement 
intervals, is responsible for the advantageous Heisenberg 
scaling when r = T [53] . 


B. Driven system with dephasing 


Real physical systems experience dephasing due to 
e.g. magnetic noise, off-resonant light scattering and 
collisions with back ground gas [s n\. To study how 
such effects affect parameter estimation, we introduce a 
transversal dephasing of the atomic coherence at a rate 7 
with the Lindblad operator in Eff- Q • The ground 

and excited states are eigenstates of Uz, and Tr {pqC[po\) 
still vanishes so that by Eq. (131 the QZE is present in 
the system. 

The system evolution may be solved analytically, and 
on resonance one finds 


P{9\9:^) 


\ -\- ^ sin UJT -\- cos LOT 

i + '^(Or + l) 


) H 7 ^ 7 
H = 7 , 

( 20 ) 


where to = ~ ■ 

Results are shown in FIG. which depicts the time 
evolution of P{g\g, H) for different values of 7 . Due to 
the dephasing, the oscillations are damped over time, and 
it is intuitively clear that it is no longer optimal to wait 
as long as possible between measurements. 


FIG. 4. (Color online) The time evolution of P{g\g,Q,) for a 
two-level system driven on resonance. Results are show for 
values of 7 ranging from no dephasing (y — 0) to critical 
damping of the Rabi oscillations (y — fl). 


Using Eq. (16 1 we determine the Fisher information 
per measurement 

Fr(il) —(qwr coswT — (7 -I- sinwr)^ 

N — 2a;^e^'>"^ -I- — 7 ^) cos 2wr -|- 2yLo sin 2lot) ’ 

( 21 ) 


valid for H 7 ^ 7 as well as in the limit H —>■ 7 . 

In the six upper panels of FIG. the Fisher informa¬ 
tion per time for different dephasing rates is depicted. 
Note the oscillatory behaviour, and that it is indeed 
peaked for finite measurement intervals. 

From Eq. (21 1 one finds, in the strong driving regime. 


Fr{0) sin^wr 

- ' 7"0 ' - 

T 1 — cos^ LOT 


for H ^ 7 . ( 22 ) 


For small r the envelope yields the same linear 

dependence as in the undamped case caused by the 
Quantum Zeno effect, while for larger r dephasing 
suppresses the information. 


The lower panel of FIG. shows in color the Fisher 
information as a function of 7 and r, and the black curve 
indicates the values of r that maximize the information 
for any given 7 . For high 7 , phase information is only 
significantly maintained during the first tt/2 of the rota¬ 
tion and quite frequent measurements are optimal. As 7 
decreases, additional tt rotations are favourable between 
measurements, and since the sensitivity to H by popu¬ 
lation measurements is highest when both populations 
are close to i, the optima occur around lot = 7r(n + \) 
where n is an integer. The optimal measurement intervals 
change in steps as the dephasing rate becomes smaller 
and smaller. 


VI. BAYESIAN INFERENCE 

In the sections above we investigated the theoretical 
optimal sensitivity based on projective measurements. 
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FIG. 5. (Color online) The Fisher information Eq. (16 1 per 
time for Rabi frequency estimation in a two-level system 
driven on resonance. In the six upper panels, results are 
shown as a function of the duration of the intervals between 
measurements and values of 7 ranging from no dephasing 
(j = 0) to critical damping of the Rabi oscillations (7 = fl). 
The color plot pictures the Fisher information as a function of 
both T and 7 . The color scale is in units relative to the maxi¬ 
mum value. The black curve in the plot tracks the value of t 
maximizing the Fisher information for any given value of 7 . 
For weak damping, the maxima occur close to cur = 7 r(n+ |), 
where n is an integer. 


We shall now address the achievement of this sensitivity 
by an explicit estimation strategy. Given the outcome 
at the time tj=jT, Bayes rule Eq. ([^ yields an update 
for the probability P{9) assigned to different values of 9 
prior to the detection event. The update rule is applied 
at each detection and thus acts as a filter multiplying all 
the probability factors assigned to the outcomes actually 
occurring during the experiment on the prior, say, uni¬ 
form probabilities assigned to the candidate values for 
the parameter 9 [15] . 

We illustrate this principle in the lower panel of FIG. 
by propagating probabilities based on three candidate 
values of the Rabi frequency used to generate the tra¬ 
jectory in the upper panel. The trajectory is dominated 
by intervals of projections to the same state in consecu¬ 
tive measurements (the QZE), consistent with a vansh- 











A 





FIG. 6 . (Color online) The top panel shows the ground state 
population resulting from a quantum trajectory simulation 
with If = flo and 7 = 5 = 0. The lower panel shows the 
evolution of the probabilities P{fl\D) for three candidate val¬ 
ues of the Rabi frequency, conditioned on the measurement 
record in the upper plot. Results are shown for flgr = 1. 


ing Rabi frequency and hence causing an increase in the 
probability for the lowest value of il. Due to the finite 
value of r, however, the system performs quantum jump 
like dynamics between the eigenstates at a rate ~ 
favoring the higher values of Ft in the filtering process. 
Eventually, after a long sequence of measurements, the 
full record is only consistent with the true value of Ft. 

FIG. shows the results of applying the same proce¬ 
dure for a wider range of candidate values of the Rabi 
frequency on a fine grid and for a longer time. A nat¬ 
ural estimate for the unknown parameter 9 is the value 
maximizing the resulting P{9\D), and with a uniform 
probability distribution prior to the measurements such a 
maximum likelihood estimate actually saturates the GRB 
Eq. ([ 6 ) in the asymptotic limit of many measurements 
m- To distinguish infinitesimally different hypotheses 
(parameter values) the probing intervals should be cho¬ 
sen to maximize the Fisher information, and with a finite 
dephasing 7 = 0.1 Ho, this happens for Dor = 4.83 (See 
FIG. [^. However, as seen in (a), the procedure is not 
able to filter away other candidate values Ft fulfilling 

\/FP -f 7 ^r = 2n7r ± ^Dq -|- (23) 


with n an integer. If D ^ 7 , these Rabi frequencies 
lead to the same ground (excited) state populations at 
the time r, and thus yield identical probabilities for the 
measurement outcomes. 

A solution to this problem is analyzed in FIG. |7b[ 
where shorter time intervals with FIqt = 0.3 are used 
to distinguish the discrete frequencies in Eq. (23l. As 
observed in the upper panel, however, we are here far 
in the Zeno regime, and even though the distribution 
centers around the true value in the lower panel, it is 
broader than the individual peaks in (a). 
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FIG. 7. (Color online) The upper panels show the ground state population during a simulation with fl = Ho, 5 = 0 and 
7 = O.lflo- The lower panels shows the evolution of the quasi continuous probability distribution P{fl\D) for the Rabi frequency, 
conditioned on the measurement records in the upper plots. Results are shown for (a) flor = 4.83 and (b) Qor = 0.3. In (c) the 
measurement is optimized (see text) with 12 measurements with a separation given by HoT" = 1-82 followed by 16 measurements 
with a separation given by fior = 4.83. The dashed, black lines represent values ofQ fulfilling the condition in Eq. (23l. Note, 
that is assumed to lie between 0 and 2.5fIo. 


An optimized strategy applies a combination of q 
measurements with sub-optimal evolution times to 
distinguish the discrete frequencies in Eq. ( [2^ , i.e., 
to select the correct branch in FIG. |7b| followed by 
N — q measurements with the optimal evolution intervals 
Toptimai) maximizing the Fisher information. From 
Eq. (231 all problematic candidates are ruled out by 


setting Ts < {2tt - ?7)/(v^^^max + 7^ + +7^), where 

f^max is the maximal candidate and ry is a small positive 
number. 

This procedure is equivalent to one used in state esti¬ 
mation m, where an asymptotically insignificant frac¬ 
tion q = N^~'^ of the measurements is applied to identify 
an estimate for the state within a neighbourhood of size 
around the true state po- The state pnir) at 
time T corresponds to a particular value of Tl, and with 
a distance from the true state given by 


L = ||pn(Ts) - Po(rs)||i, 


(24) 


where ||/x||i = Tr)^/ pip), efficient distinction of fl from 
rio by the first q measurements, dictates a lower bound 
for e, 


1 InL 

^ - 2 hTiv’ 


(25) 


We pick L to be small enough to discern values of 
the Rabi frequency that would lead to opposite measure¬ 
ment probabilities, i.e., for ^ 7 , we set in Eq. (24) 
= tt - + so that P{g\g,Tl) = 

1 — P{g\g,TtQ). Note, that L can be rewritten as twice 
the difference in ground state populations. 

Finally, we obtain the optimal parameters for N, q 
and Tg numerically subject to the total time constraint 


P — Q'^s T (A^ Q)T)ptimal- 

This scheme is implemented in FIG. The opti¬ 
mized measurement parameters are given in the figure 
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FIG. 8. (Color online) Illustration of the Quantum Zeno ef¬ 
fect for a two-level system driven on resonance subject to con¬ 
tinuous dispersive probing with strength F of the eigenstate 
populations, (a): F = 0.05II, (b): F = 0.2II, (c): F = II and 
(d): F = 312. Note, the resemblance with the plots in FIG.[^ 


caption. The measurements in the first t ~ 22f2“^ effi¬ 
ciently filter away the problematic candidates of Eq. (231, 
and the subsequent optimal measurements provide a nar¬ 
row peak around the correct Rabi frequency. Hence, the 
hybrid method proves successful. Note that while we 
presented the scheme as if the first measurements are 
needed to identify the approximate range of H first, it 
was recently demonstrated how adaptive measurements 
relying on real time Bayesian updating can be applied for 
optimized sensing with a single electron spin |25| . 


VII. CONCLUSION AND OUTLOOK 


In this article, we have presented an analysis of the 
Fisher information for the determination of parameters 
that govern the dynamics of a quantum system subject to 
projective measurements at a given frequency. We iden- 































































tified and analyzed the inhibition of precise parameter es¬ 
timation cansed by the Quantnm Zeno effect. Our results 
directly relate the Zeno inhibited evolution of a quantum 
system to a vanishing information acquired in any given 
time from the measurements causing the inhibition. Our 
exemplification through the monitored Rabi oscillations 
of a two-level system shows that, in the Zeno limit of 
frequent probing, the parameter sensitivity approaches 
zero, and that the interval of propagation between mea¬ 
surements should be made as large as possible. In the 
presence of dephasing, finite intervals are favourable for 
optimal estimation purposes, while the Zeno effect still 
excludes frequent probing. 

We emphasize that in the absence of a Zeno effect, 
where the linear term in the survival probability Eq. (11) 
is non-zero, the Fisher information per time may still van¬ 
ish in the limit r —>■ 0 if this term is independent of 9. As 
an example, consider the inclusion of spontaneous decay 
at a rate 7 spont in our two-level model. With an initial ex¬ 
cited state this yields Tr(po^[po]) = 7spont ^ 0 , and the 
Fisher information for estimating Hamiltionian parame¬ 
ters indeed vanishes, while the information for estimating 
7 spont remains finite in the limit r —>■ 0 . 

The Fisher information and the Cramer Rao bound 
quantify the optimal asymptotic sensitivity, but to reach 
unambiguous sensing at that limit, we must employ 
a strategy that distinguishes among remote candidate 
parameter values. We discussed and demonstrated a 
Bayesian estimation protocol, that spends a small frac¬ 


tion of the total probing time to establish a rough esti¬ 
mate, from which the asymptotic results are valid. The 
problem of parameter estimation from periodic probabil¬ 
ity distributions is elaborated in HZ]. 

For simplicity, and for the possibility to establish ana¬ 
lytical results, our model focused on projective (strong) 
measurements performed at discrete times, and we in¬ 
troduced the Zeno effect as the limit of very frequent 
measurements. We believe, however, that our results are 
also valid for the interpretation of different measurement 
scenarios, such as continuous, dispersive probing. In such 
experiments, a system is interrogated by a laser or mi¬ 
crowave field, which undergoes a phase shift that depends 
on the state of the system, and by adjusting the inten¬ 
sity of the probe field, both weak and strong (projective) 
measurements may be implemented on the time scale of 
the system dynamics. The probe strength T thus plays 
a role equivalent to the frequency 1 /r of projective mea¬ 
surements in the present work, cf., the similarity of the 
conditional dynamics in FIG. for continuous probing 
simulated according to Ha and FIG. 0 for the fre¬ 
quent projective measurements. 
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